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We consider a one-dimensional persisent random walk viewed as a deterministic process with 
a form of time reversal symmetry. Particle reservoirs placed at both ends of the system induce a 
density current which drives the system out of equilibrium. The phase space distribution is singular 
in the stationary state and has a cumulative form expressed in terms of generalized Takagi functions. 
The entropy production rate is computed using the coarse-graining formalism of Gaspard, Gilbert 
and Dorfman. In the continuum limit, we show that the value of the entropy production rate is 
independent of the coarse-graining and agrees with the phenomenological entropy production rate 
of irreversible thermodynamics. 
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I. INTRODUCTION 

One of the most interesting and, at the moment, controversial problems in non-equilibrium statistical mechanics is 
to understand the mircoscopic origins of the entropy production in non-equilibrium processes, especially the entropy 
production envisaged by irreversible thermodynamics for the usual hydrodynamic processes. By this one means that 
one should start from some totally microscopic description of a system in some intitial state, show that in the course 
of time the system approaches a stationary state, in some sense, and then calculate the entropy production associated 
with this process. There is no clear definition of the microscopic entropy production, however, and this in itself 
presents a problem. One choice of an entropy was provided by Gibbs with a definition based upon the full phase 
space distribution function of a classical system. However, the Gibbs entropy defined with respect to the phase space 
distribution function remains constant in time, if the time dependence of the distribution function is determined by 
the Liouville equation, as it is for conservative, Hamiltonian systems. One solution to this particular problem has 
long been discussed, the use of the Gibbs entropy as a measure of a non-equilibrium entropy requires that every 
possible trajectory in phase space (except for a set of measure zero) be followed in time with infinite precision, but if 
one relaxes this requirement and only follows trajectories to within some specified precision, i.e. "coarse grains" the 
description, then one obtains an entropy function that increases with time. Of course, an increase with time does not 
automatically imply an agreement of the entropy production with the laws of irreversible thermodynamics. 

As a consequence of these considerations, a number of issues remain to be resolved: 

1. Is a coarse grained Gibbs entropy the best candidate for a definition of a non-equilibrium entropy ? 

2. If so, how does one correctly define the coarse graining process ? To what extent are results so obtained 
independent of the coarse graining procedure ? 

3. Do any definitions of entropy production lead to the laws of irreversible thermodynamics, and if so, when and 
why ? 

There is, of course an enormous literature on all of these questions. Here we wish only to discuss some recent 
progress in answering them based upon the approach to a theory of irreversible processes through dynamical systems 
theory. Gaspard |(J has considered a microscopic model of diffusion of particles in one dimension, called a multi-baker 
model, where the dynamics is modeled by a baker's transformation taking place on a one-dimensional lattice, where 
a unit square is associated to every site. Here the one-dimensional lattice is identified as the configuration space 
of a random walker where diffusion takes place. The baker's transformation exchanges points of the unit squares 
between neighboring cells and allows for a deterministic description of the random walk, i. e. keeping the dynamics 
on the configuration space unchanged. The variables on the unit square are irrelevant to the diffusion process and, 
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in themselves, have no physical meaning, other than insuring the measure-preserving nature and reversibility of the 
dynamics. Gaspard considered a finite length chain, 1 < n < L, and defined the dynamics on the unit squares so that 
particles would be sent either to the adjacent right or left intervals depending upon their location along the expanding 
direction in the unit square. Tasaki and Gaspard [jl"if considered the case where a steady gradient in particle density 
was maintained along the chain, and were able to show that fractal-like structures formed by regions of differing 
microscopic densities, appear in the two dimensional phase space. The fractal like structures become real fractals in 
the infinite volume limit, as L — > oo. 

These fractal-like structures are strict consequences of the dynamics given the presence of a density gradient 
produced by particle reservoirs at the boundaries. Their importance for the theory of entropy production, as pointed 
out by Gaspard lies in the fact that they provide a fundamental reason for coarse graining the distribution 

function. That is, for large systems, the microscopic variations in density take place on such fine scales that no 
reasonable measurement process would or should be able to detect these variations. Gaspard |t],|8| showed that the 
steady state production of entropy in this model is in agreement with the predictions of irreversible thermodynamics 
and further, that the entropy production is independent of the coarse graining size over a wide range of possible 
coarse graining sizes. While it is not entirely clear why this procedure leads to results in agreement with irreversible 
thermodynamics, the model and procedure are sufficiently interesting and stimulating that one would hope to find 
further examples so as to gain some deeper insights into the nature of entropy production, hopefully in general, and 
certainly in this group of baker transformation like models. 

A closely related, independent approach to the problem of entropy production in non-equilibrium steady states is 
provided by Tel, Vollmer, and Breymann (TVB) in a series of papers p],|l9|,p|]2"c|] , also devoted to diffusion in multi- 
baker models. These authors considered the entropy production in measure preserving maps as well as in dissipative 
maps that do not preserve the Lebesgue measure and model systems with Gaussian thermostats. The TVB systems 
also show that a coarse grained distribution function leads to a positive entropy production in agreement with non- 
equilibrium thermodynamics. In addition to considering a wider class of models than Gaspard, they used a different 
coarse graining scheme that was not devised to expose the underlying fractal structures of the SRB measures of 
the two dimensional phase space associated with their models. They also argued that their results for the entropy 
production should be largely independent of the type of coarse graining scheme used. 

A generalization of the methods of Gaspard (?]|| and TVB [Qjlj|||j20| was proposed by the present authors in a 
recent paper Q, where it was shown that a coarse-grained form of the Gibbs entropy, which can be expressed in 
terms of the measures and volumes of the coarse graining sets partitioning the phase, leads to an entropy production 
formula similar to Gaspard's and applicable to more general volume-preserving as well as dissipative models, such as 
those considered by TVB, as well as multi-baker maps with energy flow considered by Tasaki and Gaspard ]lq] . 

We should also mention that multi-baker models and the use of coarse gaining methods for calculating their entropy 
production have been criticized by Rondoni and Cohen fllSfl . While some of their points are indisputably correct, their 
approach does not suggest a better way to proceed. So we continue this line of research in the expectation that it 
will lead to some further insights into entropy production in more realistic models, despite the shortcomings of the 
simplified models we treat here. 

In this paper, we extend the ideas mentioned above, particularly those of Gaspard Jt],|| and Gilbert and Dorfman ||, 
to a somewhat more complex model than those considered above. We consider a deterministic version of a persistent 
random walk in one dimension. The persistent random walk (PRW) is similar to the usual random walk , but in this 
case the moving particle has both a position and velocity as it moves along a one dimensional lattice. At each lattice 
site the particle encounters a scatterer which, with probability p allows the particle to continue in the direction of its 
velocity and with probability q = 1 — p reverses the direction of the velocity. This model is a limiting version of a 
Lorentz lattice gas described in detail by van Velzen and Ernst jl8j , where scatterers are distributed at random along 
the lattice sites with some overall density per site. In the persistent random walk, the site density is unity. Although 
this is clearly a random process, it can be turned into a deterministic one by a method described by Dorfman, Ernst, 
and Jacobs based upon the baker's transformation, whereby new variables are added to the system such that the 
dynamics in terms of these new variables allows us to replace the stochastic scattering mechanism by a deterministic 
one. A similar mechanism has also been considered by Goldstein, Lanford and Lebowitz in a different context p0[ . 
Further, we can place many particles on the lattice as long as they do not interact with each other and if we take 
each individual scattering event to be independent of any others taking place at the same instant of time. 

Here we will consider the PRW in one dimension as a model for diffusion and entropy production, and we will describe 
the entropy production in terms of the additional phase space variables needed to make the system deterministic. 
That an irreversible entropy production is to be associated with this process follows from the fact that for systems 
with periodic boundary condiditons, any intitial variations in the probability of finding a particle at given positions 
will eventually vanish and the probability will become uniform with time. Thus information is lost about the position 
of the particle and the entropy of the system is thereby increased. It is of some interest to see how this loss of 
information is reflected in the phase space distribution function, and how the entropy so produced is related to the 
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entropy which is the subject of irreversible thermodynamics. We will see that the analysis of the PRW model has 
some features that one hopes would be more general, such as the clear independence of the entropy production on the 
sizes of the coarse graining regions, provided these regions are not too small, and the coincidence of the production 
of the coarse grained Gibbs entropy with the entropy production of irreversible thermodynamics, previously noted by 
Gaspard, TVB and the present authors in multi-baker maps. 

The organization of this paper is as follows. Section || gives a simple description of the phenomenological approach 



to entropy production in a simple one-dimensional diffusive system. In Sec. [II, we describe in more detail the 
PRW process. In Sec. [V we describe a two-dimensional, measure preserving map that reproduces the PRW on a 
macroscopic scale and discuss its time-reversal properties. We consider a system of non-interacting particles with a 
steady density gradient, produced by particle reservoirs at each end of the system. We then obtain the steady state 
invariant measure using techniques due to Gaspard and Tasaki |fl4||. In Sec. |v| we develop a symbolic dynamics for 
the diffusion process, as it is refected in phase space, and, in Sec - [VJ, we apply this symbolic dynamics to compute 
the rate of entropy production. We conclude with remarks and a discussion of open questions in Sec. VII. 



II. PHENOMENOLOGICAL APPROACH 



Before turning to the description of the persistent random walk in Sec. [II, we briefly discuss the phenomenological 
approach to entropy production in a one-dimensional system. 

Let us consider a stochastic diffusive process that is driven away from equilibrium. More specifically, we have 
in mind a system where a gradient of density is imposed by appropriate boundary conditions. Let the system be 
essentially one-dimensional by imposing translational invariance in the other spatial directions. The relevant spatial 
direction is denoted by x, < x < L + 1. At the boundaries, x = and x = L + 1, we put particle reservoirs with 
respective particle densities w_ and w + . The reservoirs keep on feeding the system with new particles and those 
particles that exit never come back. 

The probability density of a tracer particle is Wt{x) and obeys the mass conservation law 

dw t (x) 



dt 



+ V-j t (aO = 0, (1) 



where jt{x) is the associated current. 

To follow the phenomenological approach of non-equilibrium thermodynamics, we have to supplement the mass 
conservation law with a linear law that relates the current jt(x) to the gradient of the density Vwj(x), 

j t (x) = -DVw t (x), (2) 

where D is the diffusion coefficient associated with that process. This is known as Fick's law j|. 

Combining Eqs. (|l|-|^) and assuming that the diffusion coefficient is constant, we obtain the Fokker-Planck equation 
for diffusion, 

dwt(x) 2 

dt = D V w t (x). (3) 

The stationary solution of this equation is found by imposing the boundary conditions Wt(0) = u>_ and w t (L+l) = w + , 

x 

w(x) = W- + (w+ — w-) — . (4) 

L + 1 

The connection with the second law of thermodynamics is made by considering the entropy whose local density is 
defined as 

a t (x) = -{logMz)] - 1}, (5) 

and S(t) — Jq +1 dxwt(x)s t (x) is the macroscopic entropy. Taking the derivative of the integrand with respect to 
time, we find successively 

dw t (x)s t {x) dw t (x) . . 

— m — = -^r log[wt{x)] 

= [V- j t (x)]log[w t (x)] 

= V-{j t (x)log[w t (x)}} + ^^, (6) 
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where we used Eq. (|I|) in the second line and Eq. (^J) in the third one. Equation (g) has the form of a local entropy 
balance || 

dw t (x)s t (x) _ not 



dt 



= -V-Jl ot ( Xl t)+a t (x), (7) 



where Jl ot (x,t) is the total entropy flow and Ot{x) > is the entropy source term. Equation (|7|) can be rewritten in 
a slightly different form 

w t 0r)^P =-V- J s (x,t)+a t (x), (8) 
where the entropy flux is the difference 

J s (x,t) = Jl 0t (x,t)-w t (x)s t (x)j t (x), (9) 

between the total entropy flux and a convective term. 

Using Eq. (^|) again, we find that the local rate of entropy production <7t(x) in Eq. (|^) can be rewritten in terms of 
the density Wt(x) only : 

M-D&S&f. ( 10 ) 

Wt(X) 

In the stationary state, Eq. (Q), the local rate of entropy production becomes 

a(x)=D {w+ -™- )2 . (11) 



III. PERSISTENT RANDOM WALK - DEFINITIONS 



We consider the following process : a particle on a one-dimensional lattice moves from site to site with a given 
velocity ±1 and, at each time step, is scattered forward with probability p and backwards with probability q, such 
that p + q = 1. Thus a particle located at site n S Z with velocity v = ±1 will move to n ± 1 conserving or reversing 
its velocity with the probabilities : 



(n,v) 



(n + v, v) 
(n — v, ~v) 



with probability p, 
with probability q. 



(12) 



The diffusion coefficient for this process is flqjlqjl 



D 



P_ 

2q 



(13) 



To construct a deterministic model of this process, we start with a procedure similar to the one used by Dorfman, 
Ernst and Jacobs and consider the velocity of the particle only and associate to it a unit interval. Let us divide 
the interval into two halves, corresponding to the two possible values of the velocity. For definiteness, let us assign 
the first half to v = — 1 and the second to v = +1. According to Eq. (|l2|), a particle has a probability p of keeping 
its velocity unchanged and q to reverse it. Therefore we have to further subdivide both halves into two parts of 
respective lengths p/2 and q/2, the first of which must be mapped onto the same half and the second onto the other 
half. Each of these branches must be linear and onto in order to model the independent process of velocity flips. This 
is illustrated in Fig. (|l]) for the case q — 1/3. Therefore a one-dimensional map whose dynamics mimics the random 
sequence of the velocities is given by 



<f) q (x) 



( x , 1 

p ' 
x - H 



< x < |, 



§<*<!-§, 



(14) 



< x < 1. 
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The generalization of &q(x) to a multi-baker map is straightforward 




1,£ 



1 



2>« W-2 



1 



n — 1 , tpt 1 — 



g - g/2 
1 x ~ 1 + 



g/2 



X-1 + q/2 




0<x<§, 0<y< ^, 
0<z<§, ^<y<l, 



< a; < 



< x < 1 - 



(15) 



3 y 



%<x<l,0<y<%, 
%<x<l, \ <y < 1. 



This map is shown in Fig(||). We note that if one restricts this map to a single unit square, by ignoring the changes in 
the site index, n, this map is reversing under the time-reversal operators of the usual baker map T(x, y) = (l — y,l — x) 
and S(x,y) = (y,x), i. e. 



ToMV™oT(x,y) = MP™- 1 (x,y), 
SoM»™oS(x,y) = M q >™- 1 (x,y),. 



(16) 
(17) 



However, the full multi-baker map, with the changes in the site index taken into account, reversing under either of 
these two time-reversal operators. Indeed, one finds that, depending on x and y, the composition 



T o Mf w o T o Mf w (n, x,y) = 




,x,y 



(18) 



and similarly for S. This is displayed in Figs. (^Q). It is interesting to note that we can combine the action of S and 
T to form a hybrid operator U that is reversing. Indeed, as seen in Fig. (§), S and T act on the square in such a 
way that the images of the elements of the natural partition of the square under these two operators do not overlap. 
Therefore U, defined by 



' T(x,y), Q<x<\, §<y<5, 

i-f <y<i, 

U(x,y)={ \<x<\, 0<?/<|, 

l <y<i-i 



(19) 



, S(x, y) otherwise, 



is a reversing operator for M 9 prw . However, contrary to T and S, U is not an involution, i.e. U o U{x,y) ^ (x,y), and 
therefore U is not necessarily a reversing symmetry of any power of M ? prw . Actually U o U is an involution, 



U 4 (x, y) = U oU oU o U(x, y) = (a:, y). 



(20) 



This implies that the conjugation of U with odd powers of M ^ rw yields its inverse, while the conjugation with even 
powers leaves the map unchanged. In other words, U is a reversing symmetry of AfP rw ™ for n odd but a symmetry 
for n even : 



U o A/?™" o U = 



MP™"' 
MP™", 



n odd, 
n even, 



(21) 



This property of a reversing operator can be contrasted with the existence of two reversing symmetries for the maps 
on the unit square, which itself implies the existence of a non-trivial symmetry, i.e. the composition of T and S. Those 
properties are well established [|ll|, while the property described by Eq. flgj] ) is new, to our knowledge. 

The weakened reversibility of the map MP™, Eq. (|2l]), still allows the map to have a properly behaved dynamical 
entropy production as we will see in Sec. [Vl|. 



5 



IV. STATIONARY STATE UNDER FLUX BOUNDARY CONDITIONS 



As with the open multi-baker map (IJJ^,P, we impose flux boundary conditions. We consider a chain of L sites 
and study the distribution of an infinite number of copies of identical systems imposing that, in average, 1 particle is 
present at the left-end of the chain, regardless of its velocity and L + 2 particles at the right-end. 

The stationary measure is found by considering the following cumulative functions : 



&-\n,x,y) 



dx' 



v 1 

dyp(n,x,y), 0<x<-, 



G i+) {n,x,y) = dx' I dy'p(n,x',y'), 

Jl/2 



< X < 1, 



(22) 
(23) 



where p(n, x, y) is the corresponding density function. Because the map is piecewise uniformly expanding along the 
x direction, the invariant measure is uniform along the x-direction and we therefore have 



G [ -\n,x,y) = 2xg { - ) (n,y), G {+) {n,x,y) = (2x - l)<? (+) (n, y), 
where the functions and are solutions of the following set of equations : 

' qg& (n±l,f + ^) -qg<& (n ± 1, , < y < §, 

qg& (n ±1,1)- qg^ (n ± 1, J) + P9 iT) ( n ± !> ~ 



(24) 



< 



y<i-H, 



qg& (n ± 1, 1) - qg {±) (n ± 1, ^) + P9 iT) (n ± 1, 1) 
■ w (±)f n±1> l^i+l), l-j<„<l, 



9 (T) {n,y) = < 



q 1 2 j 
With flux boundary conditions, 

5 (-)(0, y) + g<-+\0, y) = y, g^(L + 1, y) + gM(L + l,y) = (L + 2)y, 
the solutions of ( |25| ) for y = 1 are readily found to be 

g^ >(n, 1) = — - — =F — . 



(25) 



So that, for arbitrary y, the solution takes the form 

9 {±) (n,y) 



n+1 1 
— ^— T -r 



±^(1/), 



(26) 



(27) 



(28) 



where the class of functions Tn is a one-parameter family, i.e. implicitly dependent on the scattering parameter q, 
of generalized incomplete Takagi functions, 



P „ T (T) /J , 1 



l i 2 -y 

2 



(±) / y-g/2 



0<y<§, 
^ s y < i - 2 , 
1 - f < y < l, 



(29) 



where 1 < n < L and the boundary conditions on the functions T„ (j/) are 

T o (+) (y) = 0, T&G,) =0. 



(30) 



In Fig. ([j|), we display the functions T„ ' for the sites n = 1,3, 5 on a chain of L = 100 sites and the value of the 
scattering parameter q = 1/3. 

We now want to simplify the argument and will make the assumption in the sequel that we are far enough from the 
boundaries so that we can ignore the finite size effects and replace the incomplete functions by their common 
limit value 



G 



P 7/ 

T ■ \ J I/ 2 " 



l<y< 



2 
2' 



(31) 



l-f <y<i, 



where we have made explicit the parametric dependence of the generalized Takagi functions T q . Figure (0) shows a 
recursive computation of T q (y), for q = 1/3. This example is very similar to the Takagi function |l3| 



T(y) 




o<y< 2, 



^T(2y-1), 



(32) 



that appears in the example of the open random walk with flux boundary conditions discussed by Gaspard and co- 
workers 14,^,^]. In particular, for q — 1/2, the case of a symmetric persistent random walk, up to a factor 2, T q is 
identical to T on the first half of the interval and opposite on the second half. Figure (||) shows T q (y) for different 
values of q ranging from 0.1 to 0.9. 



V. SYMBOLIC DYNAMICS 



In order to compute the entropy production, it is convenient to introduce a symbolic dynamics. Each half of the 
unit cell is partitioned into sets that correspond to particles that were back- or forward-scattered at the preceding 
time step. This defines the 0-partition of a cell, 



A-Sr(-) r (+) r(-) r (+h 

— l 1 ' 1 ' 1 1 ' 1 1 J i 



where 



and 



r (±) nr 4 , « = o,i, 



(33) 
(34) 

(35) 
(36) 

l}. (37) 

(38) 

An (I, fe)-partition is the collection of cylinder sets 

rfcl..,^ = r (±) n [nt_ ; Mf wi (r Wi )] , (39) 

where € {0, 1}, i — —I, . . . , k. The measure of such a set, fj,n(u>-i, . . . , Ufc), can be written, for / = 0, as 

A 5 W(w , . . ■ , u> k ) = 9 {±) (n, y{uj , . . . , u; fe + 1)) - ,g (±) (n, y(« 0l . . . , a;*)) , (40) 
where wq, ■ ■ ■ , wje-ijCOk + 1 is defined by 



r<-) = < 




0<.,<i 




r (+) = . 








r„H 


fall) 


< t/ < - or 1 

~ 2 2 


r H 




g 

- < y < 1- 


-§}• 



W0, ■ ■ ■ , Wfc-1, Wfc + 1 = 

with the further convention that y(l, . . . , 1, 1 + 1) = 1. 



w , 



,w fe _i,l, w fe = 0, 
,w fe _i + l,0 w fe = 1, 



(41) 



By looking at the y-components of M'™ in Eq. fll5[), we find that y(u>o, . . . , w^) can be computed through the 
recursion relation 
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y(wo,...,wjk) = (v{uo)y{ui,---,Uk)-{-Y ^) modi, (42) 

where we set 

"M = ( ? ' "° = ?' (43) 

I P, UJQ = 1. 



Therefore the functions T 9 in Eq. (31) can be defined directly in terms of the symbolic sequences, 

J £ (y(ui,...,u k ) - «) -gr^wi,...,^), ^o = 0, 
r,(wo,...,w*)= < p V / i\ (44) 

[ _ 2 U^ 1 ' ■••i^feJ - 2 ] +PJg( w l> ■ ■ ■ ,Ukh w o = 1- 

With the help of Eq. ( p8| ) , we can then rewrite Eq. (|4(]) as 

A ff f ± )(w ,...,Wfc) = (^4~^ T 4) A ^K-.-,^)±AT ? K...,^), (45) 

with obvious notations for Ay and AT 9 . We note here that our assumption that the cell n is sufficiently far away 
from the boundaries means in terms of the size of the cylinder sets that k must be strictly less than the distance to 
the closest boundary. In that case, we can substitute without loss of generality Tn by T q in Eq. (|28|). 
For later purposes, we note that, from Eqs. (f42| , Q), we have 

Ay(uj Q , . . . ,oj k ) = u((J )Ay(ui, ...,uJk) = JjK^i), (46) 

and 



Krr , \ J fAy(wi, . . . ,u k ) ~ q A T q (wi, . . . ,u k ), w = 0, 

I -^Ay(wi, . . . , u k ) + pAT, wi, . . . ,w k ), uj q = 1. 



(47) 



VI. ^ENTROPY AND ENTROPY PRODUCTION RATE 



Following 



we write the entropy of a (I, fc)-partition as 



and 



?g } («) 



H { n ) {U-U---;U k -l) 



tt)_(,...,a) fe _i 



log 



Z/ (±) (w^, . . .,U) k -i) 



(48) 



(49) 



Because the invariant measure is uniform along the expanding direction, the ratio of jt/i and i/W obeys the 
identity 



(^-1. ■ ■ ■ ,^fc-l) = M^H^O: ■ • ■ L ^fc-l) 

^ (±) (w_;, . . . ^ (±) (w , . . . 

Therefore the entropy is extensive with respect to the x-direction and we have 

Henceforth, we will drop the I dependence and will simply consider the fc-entropies, which we write 



Sf ) (») = -E A ^ ±) K) 



log 



.A<?i ±} (^) 



^ (±) (^ fc ) 



(50) 



(51) 



(52) 
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where we used the compact notation uj k = ujq, . . . ,ujk~i- Here 1/^*0 (wj.) is the volume of the corresponding cylinder 



set, 



fc— l 



(53) 



The entropy production rate follows by a straightforward generalization of the formalism detailed in M to Eq. (B2 



A i S k (n)=A i si-\n)+A i S l 



?(+)/ 



and 



a4 ±} H= E A 5 W(^ fe+1 )log 



(54) 



(55) 



Assuming that the stationary state is dominated by the linear part, we can compute the fc-entropy production rate 



by expanding the expressions for AiS^\n) in Eq. ( |55| ) in powers of 

AT n 



( n + 1 + 1 \A ' 



(56) 



which we further expand in powers of — ^ ^ . To first order, Eq. ( |54| ) becomes 

AiS k (n) -- 



2_ v [AT q {u k+l ) - u{Lo k )AT q {u k )] 2 
+ 1 ^ 



Ay(w fc+1 ) 

— fe+i 

Making use of Eqs. (46-47]), we readily see that this expression is independent of /c, i.e. 

[AT q {ui k+l ) ~ v{Lo k )AT q {u k )] 2 



(57) 



E 



A y(u k+ i) 

[AT q {u k ) - l /(a; fc _ 1 )AT 9 (o; fc _ 1 )] 2 



E 

[^=0,1 wv ; 

P_ 
4g • 



[AT,( W )] 2 



Therefore, 



A i S' fe (n) 



P 1 



(58) 



(59) 



The derivation of the next order term proceeds as follows. For the sake of simplifying the notations, we will write 
Eq. (p7|) as 



fl C±)= s (±)(„,i) 
We start from Eq. (55) and substitute Eq. (|||) for Ag^\u k+1 ) : 



(60) 



A4 ±) N = E A ^ ±) (^ fc +i)iog- 



^MA^K) 

E [.9i ±)A yK+i)±Ar 9 K +1 )" 

3 tt)Ay(o; fc+1 )±AT ? (a; fc+1 ) 



log 



g^Ay(u Lk+1 )±u(u> k )AT q (u Lk ) 



(61) 
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where, in the last line, we used 



v(u k )Ay(uj_ k ) = Ay(Lu k+1 ). 



Factoring g^ Ay(uj_ k+1 ) in the logarithms and expanding up to fourth order, Eq.(|6l|) becomes 



[g^ ) Ay{u k+l )±AT q (u k+l ) 
AT q (u k+l ) 



log 



1± 



g^Ay(ui k+1 ) 



log 



]T [g^Ay{u k+l )±AT q {u k+l ) 



1 ± v{u) k )AT q (w k ) 
g^Ay{ui k+1 ) 

± AT q (uj_ k+1 ) 
s£ t) Ay(w Jt+1 ) 



=F 



V {u k )AT q {ui k ) AT,K :+1 ) 2 
fl^Ayfe^x) 2 U±)Ay(o; fc+1 ) 



v{Lo k ) 2 AT q {u k ) 2 ± Ar,(w fc+1 ) 3 



AT g fe, 



■k+l) 



j^Ay(u; k+1 ) 



=F 



J^Al/K+i) 



S^AyK+x) 



which, keeping terms up to 0(1/ (n + l) 3 ), takes the form 



J2 { ± AT g (^ +1 ) T K^)AT g (o, fc ) + Ag(^+i) 2 

JT. [ 2 9n >Ay(u k+1 ) 



v(w k )AT q (ui k+1 )AT q (iJ k ) v{u k ) 2 AT q {u k f 



=F- 



gir J Ay(^ fe+1 ) 

AT ? K +1 ) 3 



± 



2 5 ( ± )Ay(^ +1 ) 
K) 2 AT 9 (^ +1 )AT 9 (c. fe ) 2 



fc+i/ 



si^Aj/K+x) 



AT g (w fc+1 



gi^Ay(u k+1 )} 2 12 U^Ay^+J 



^fAT^^AT^) 3 v {u k fAT q {ui k ) 4 



UJ 



k+l, 



■k+l, 



Therefore, the fc-entropy production rate, Eq.(B4[), reads 



AiS k (n) = 



1 



1 



1 1 



(-) 



E 



[AT q (uj_ k+1 ) - v(uJk)AT q (w k )]' 



Ay(w fe 



■k+l J 



,(+) 2 fi«(-) 2 



E 



-AT q {ui k+1 f + 3^K) 2 AT,(^ fc+1 )Ar 9 (^ fc ) 2 - 2u{u; k fAT q {u k f 

Ay(^ +1 ) 2 
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E 



AT q (oj k+1 ) 4 - 4u(uj k ) 3 AT q (uj k+ i)AT q (uj_ k ) 3 + 3v(u; k ) 4 AT q (u k Y 

A v(^k+iY 



Assuming that n is large, we can expand the ratios involving in powers of l/(n + 1). To third order, we get 

112 1 



2gM 2g^ n + 1 2g 2 (n+l) 3 ' 
1 1 4 



1 1 4 



M +) 12^ 3(n + l) 



3 ' 



(65) 

(66) 
(67) 
(68) 



Therefore, to leading order in l/(n + 1), we retrieve Eq. (p7|), which, along with Eq. (pq), gives Eq.(|59]). 
To work out the next order terms, we first note, with the help of Eq.(J62), that expressions involving both AT g (w fc+1 ) 
and AT q (w_ k ) can be transformed into expressions involving only AT q (u k ) by summing over uj k . For instance, 



v(u k ) 2 AT q (u k+1 )AT q (iJ k ) 2 _ AT g (^ fc+1 )AT,(u, fc ) 5 



Ay(^ +1 ) 2 



(69) 



Now, by definition, 



E AT ^k+i) = E [ T Mk,"k + 1) - 



Therefore, 



E 



: ^(w^i.Wfe-i + 1) - T q (uj_ k ), 
AT q {ui k ). 



/K) 2 AT g (^ fc+1 )AT g (o; fe ) 2 _ AT g (u; fe ) 3 



and, similarly, 



E 



^K) 3 AT 9 (^ fe+1 )AT,(a; fc ) 3 _ AT,( a ) 4 



A2/(^ fc+1 ) d 



A v(^ k Y 



(70) 



(71) 



(72) 



Combining Eqs. ( p6[ - p7| , |69[]7^ ), we can rewrite the last two lines of Eq. ( 55 ) as 



v Ar,(^ fc+1 ) 3 _ v AT^f 



3q(n+ 1) 
4 

3(n+l) 3 



1 — k 

v Ar g (^ fc+1 ) 4 _ v AT q (u k ) 4 



(73) 



We can work out these expressions by substituting Eq.([47|) for AT q (uj_ k+1 ) and summing over ujq. Before doing so, 
we will need the following 



E A y(^+i) = i> 

E AT 9 K +1 ) = $>T 9 H = 0, 



(74) 
(75) 
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E 



AT q (uj_ k+1 ) 



-Ay(u k ) - q AT q ( : 



P 



Ay{u k ) -pAT q {u k ) 



pAy{ui k ) 



hi 1 * 1 

4 \q P 
(k + l) P 



E 



AT q {ui k f 

Ay(ui k ) 



■Aq 



For the cubic term, we have 



E 



AT,(^ +1 ) 3 
Ay(w fe+1 ) 2 



E 



l3 



P Ay(ui k ) 
~ 3 '_L_A 

2 2 

,q p 



= ?-l---)+(p-q)Y^^- 



_p 

* l—l 



2 E^-«) < = ^3(p-«) [i-(p-«)' 



and, for the quartic term, 



E 



AT q (uJ k+1 ) 4 
Az/K +1 ) 3 



E 



[|Ay(w fe ) - qAT g (u k ) 



fAy{ui k y 



^Ay^) - pAT q {u k ) 



P d Ay(^ fe ) J 



16 



1 1 

7 + 7 



k 



' 8q 2 Aq 



(p-q) [l-ip-q)^ 1 ] 



\- AT g (^) 4 
Ay(^ fe ) 3 ' 



(76) 



(77) 



(78) 



Equation ( |65| ) combined with Eqs. |6(j|-|68, 73, 7^-[78|) yield the third order correction to the entropy production 
rate : 



AiSk(n) = 



P 1 
2q n + 1 
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3 " 


n) 



+ 12 + 3q 2 + 2q 2k + iP q) 6g 3 



(n + 1) J 



1 



(n + l) £ 



(79) 



We note that, in the symmetric case, p = q = 1/2, the coefficient in front of the third order term is 7/12 + k/2. The 
term linear in k is identical to the case of the random walk |7],||], but the first term is different. For q ^ p, we find a 
new term, proportional to (p — q) k , which decays exponentially. Thus, for large k, Eq. ( f79| ) is at least qualitatively 
similar to the case of the multi-baker map with a linear divergence in k that is third order in the small parameter 
l/(n + l). 
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VII. DISCUSSION 



We have shown that the entropy production formalism of Gaspard, Gilbert and Dorfman applies successfully to a 
persistent random walk driven away from equilibrium by a density current. The leading order term in the expansion 
in inverse powers of the spatial coordinate is in exact agreement with the phenomenological entropy production, Eq. 
(|i~i"|). It is particularly important to note that this term is independent of the coarse-graining parameter, k. Our 
result is therefore similar to Gaspard's for the simpler case of a random walk , and we also find 

lim lira lim f n A,5 fc = D, (80) 

fc^oo |v Al „|/ Mll ^o i — 00 (V/i n ) 

where we wrote (i n = g(n, 1) and V denotes the density gradient with respect to the lattice coordinate. The limit on 
the resolution parameter is taken last, which is to say that the resolution dependent terms are accounted for by finite 
size effects and have thus no counterpart in thermodynamics, where one assumes that systems have infinitely many 
degrees of freedom. 

It is interesting to note that the form of the first order contribution to the fc-entropy production rate, Eq. (p^), is 
rather universal. Indeed similar expressions arise in other baker map models |Q,^|,^5). It would therefore be interesting 
to understand in a more general setting the connection between the form of the generalized Takagi functions and the 
expression of the entropy production involving the diffusion coefficient. This question seems to be at the heart of the 
agreement between the dynamical and phenomenological approaches to entropy production and still needs further 
explanation. 
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FIG. 1. The map <j> q defined by Eq. (|14j) mimics the probability rules, Eq. for the velocity vector v. 
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FIG. 4. The composition of Aif rw with 5, similar to Eq. (||). 
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FIG. 5. The action of 5 and T on the elements of the partition induced by M? rw . 
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